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A STUDY  OF  FLUSHING  IN  THE  DELAY/ ARE  MODEL 


Introduction 

Late  in  1951  the  Chesapeake  Bay  Institute  was  asked  to  assist  the 
Oceanographic  Division  of  the  Hydrographic  Office  in  a study  of  the  flush- 
ing problem,  using  the  Delaware  River  and  Bay  Model  located  at  the  U.  S. 
Army  Engineers  Waterways  Experiment  Station,  Vicksburg,  Mississippi. 
The  first  test  made  showed  that  the  dye  used  as  a tracer  contaminant  was 
adsorbed  in  large  quantities  by  the  concrete  bottom  of  the  model,  the 
glass  sample  bottles,  and  even  by  the  dust  which  settled  from  the  air  into 
the  model.  Further  tests  were  delayed  until  February  1952  while  the  model 
was  painted  with  adsorption  inhibiting  plastic  paint. 

Eight  separate  tests  were  made  on  the  model  by  introducing  a known 
volume  of  dye  and  tracing  the  changes  in  the  distribution  of  the  dye  by 
photometric  analysis  of  samples.  Considerable  reduction  of  data  was 
necessary  to  compensate  for  the  adsorption  and  the  body  of  this  report 
will  deal  only  with  the  reduced  data. 

Several  appendices  have  been  included.  Appendix  I deals  with  the 
sampling  gear,  the  method  of  sampling,  and  the  methods  of  field  analysis. 
Tables  of  observed  dye  concentrations,  by  station  and  tidal  cycle  are 
given  in  Appendix  II.  Appendix  III  gives  the  reduction  methods  used  to 
compensate  for  the  adsorption  and  tables  of  corrected  concentrations  by 
station  and  tidal  cycle.  Appendix  IV  contains  auxiliary  data  on  salinity, 
high  and  low  water  profiles,  and  fresh  water  inflow  as  well  as  tables  of 
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cross'-sectional  areas,  inter  sectional  volumes,  and  accumulated  volumes 
for  each  5000  foot  station. 

Description  of  Flushing  Tests 

The  Delaware  Model 

The  region  represented  by  the  Delaware  model  includes  all  of  the 
Delaware  Bay  and  the  tidal  section  of  the  Delaware  River,  from  the  Capes 
to  Trenton.  The  model  is  constructed  of  concrete  on  a horizontal  scale 
of  1 to  1000  and  a vertical  scale  of  1 to  100.  A chart  of  the  model  is 
given  in  Figure  1. 

Information  about  the  hydraulic  adjustment  and  verification,  and 
the  salinity  verification  are  presented  in  various  reports  of  the  Water- 
ways Experiment  Station  (Corps  of  Engineers,  U.  S.  Army  1951(a), 
1951(b),  and  1952).  These  reports  were  sufficiently  conclusive  to 
suggest  that  valuable  information  or>  flushing  in  the  Delaware  River  and 
Bay  could  be  obtained  from  the  model. 

The  Flushing  Te3ts 

In  the  following  discussion  the  tests  are  grouped  according  to  the 
point  of  release.  In  Tests  1=4.,  1,  and  4 the  dye  was  introduced  at  Station 
52.5  which  is  located  a short  distance  below  Philadelphia.  Station  292.5 
below  Artificial  Island,  was  used  as  the  release  point  in  Tests  2 and  3. 

The  dye  was  released  from  Station  111  located  near  Chester,  for  Tests 
5,  6,  and  7. 

While  some  treatment  is  presented  of  the  lateral  and  vertical  distri- 


butions, the  main  discussion  here  will  center  on  the  variation  in  concentration 
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along  the  axis  of  the  estuary.  In  each  case  the  data  represent  conditions 
at  the  time  of  high  water  slack. 

The  manner  of  introducing  the  dye  is  discussed  in  detail  in  Appendix 
I.  Briefly,  the  dye  was  introduced  at  high  water  slack  in  a cylinder  con- 
structed so  that  the  sides  could  be  freed  from  the  bottom  and  pulled  out 
of  the  water,  leaving  the  dye  suspended  momentarily  in  the  water  in  a 
cylindrical  9hape.  Then  as  the  ebb  tidal  movement  started,  the  dye  was 
rapidly  spj  ead  in  both  a lateral  and  a longitudinal  direction.  The  diameter 
of  the  cylindrical  container  was  0.74  ft,  equivalent  to  740  ft  in  the  proto- 
type. The  height  of  dye  in  the  container  was  adjusted  to  be  even  with 
the  surface  of  the  water  in  the  model  at  the  time  of  release. 

Test  1 -A,  1 , and  4 - Release  Point  52.5 

Test  1-A,  which  was  the  trial  run  made  in  December  1951,  and  Test 
1,  which  was  the  first  of  the  series  of  runs  made  in  February  1952,  were 
duplicate  tests,  the  dye  being  released  at  Station  52.5  under  conditions 
of  mean  tide  and  mean  river  flow  (12,350  sec  ft  at  Trenton).  The  dye 
concentration  at  release  was  1000  ppm,  and  even  though  at  this  location 
in  the  tidal  river  the  water  is  practically  fresh,  no  density  effects 
were  observed. 

Upon  release  the  dye  cylinder  moved  seaward  with  the  ebbing  tide 
while  mixing  laterally  and  longitudinally.  After  one  tidal  cycle  the  dye 
had  spread  longitudinally  over  nearly  one  tidal  excursion  assuming  a 
nearly  normal  distribution,  and  after  two  tidal  cycles  it  had  a nearly 


uniform  distribution  laterally. 

It  was  observed  that  in  pockets  in  the  shore  line  small  eddies 
would  sometimes  trap  water  containing  a high  concentration  of  dye. 

After  the  high  concentration  of  dye  had  passed,  the  entrapped  dye  would 
be  slowly  fed  back  into  the  main  stream.  This  phenomenon  was  observed 
both  on  the  ebb  and  on  the  flood  tide  and  probably  materially  assists  the 
spread  of  the  highly  concentrated  contaminant.  Figure  2 is  a schematic 
representation  of  the  process. 

Figures  3a  and  3b  show  the  longitudinal  distributions  of  the  con- 
taminant concentration  from  1 to  5 8 tidal  cycles  after  the  dye  was  released 
for  Test  1-A.  (Test  1 was  substantially  the  same  as  Test  1-A.)  The 
solid  vertical  column  shown  at  Station  52.  5 represents  the  initial  dye 
volume.  After  one  tidal  cycle  the  contaminant  had  spread  over  nearly 
one  tidal  excursion  and  the  distribution  approximated  the  normal  curve. 
With  increasing  time  the  spread  increased  while  the  peak  concentration 
moved  seaward. 

The  rate  of  decrease  of  the  peak  concentration  seems  to  be 
exponential  and  Figure  4 shows  a semilog  plot  of  peak  concentration 
against  time  in  tidal  cycles  for  both  Test  1-A  and  Test  1.  Two  expon- 
ential decay  rates  are  evident.  For  the  first  two  tidal  cycles  the  rate  of 
decrease  is  extremely  rapid,  falling  from  1000  ppm  for  time  zero  to  less 
than  6 ppm  after  two  tidal  cycles.  For  times  greater  than  two  tidal  cycles, 
decay  progresses  at  a much  slower  rate.  It  would  seem,  then,  that  the 
process  governing  the  dispersion  of  the  contaminant  during  the  first  two 
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tidal  cycles  are  fundamentally  different  from  those  acting  after  the 
second  tidal  cycle. 

Figure  5 shows  the  movement  of  the  peak  concentration  with  time 
in  tidal  cycles  by  the  solid  line.  The  dashed  line  is  the  calculated  net 
downstream  water  movement  based  on  river  inflow.  It  is  to  be  noted  that 
the  peak  concentration  moves  downstream  more  slowly  than  the  net  down- 
stream water  movement.  This  was  the  case  in  every  test. 

If  instead  of  following  the  movement  of  the  peak  concentration,  we 
consider  the  concentration  of  dye  at  a fixed  point,  we  get  a different 
picture.  Figure  6 shows  the  change  in  concentration  with  time  for  sta- 
tions above  the  release  point  while  Figure  7 does  the  same  thing  for 
stations  below  the  release  point.  For  stations  well  above  the  release 
point  (e.g.  20  and  30)  an  increase  in  concentration  occurs  during  the 
first  few  tidal  cycles  and  is  followed  by  an  approximately  exponential 
decay  after  the  local  peak  has  been  passed.  For  upstream  stations 
nearer  the  release  point  (e.g.  40  and  50)  the  initial  rise  occurs  during 
the  first  tidal  cycle.  The  decay  is  very  rapid  during  the  first  few  tidal 
cycles  becoming  approximately  exponential  as  time  passes.  For  stations 
downstream  from  the  release  point  the  decay  becomes  approximately 
exponential  with  the  passage  of  time  but  the  initial  increase  prior  to  the 
passage  of  the  peak  concentration  is  not  a simple  exponential  function. 

The  rate  of  increase  in  local  concentration  as  the  peak  concen- 
tration approaches  becomes  progressively  less  for  stations  further  from 
the  release  point.  This  can  be  accounted  for  by  the  fact  that  as  the  peak 
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concentration  approaches  a given  station  the  size  of  the  peak  is  decreasing. 
So  long  as  the  rate  of  approach  of  the  peak  concentration  is  faster  than  the 
rate  of  decay  of  the  peak,  the  local  concentration  will  increase  with  time. 
For  stations  well  below  the  release  point,  the  decrease  in  peak  concen- 
tration and  the  increase  in  spread  may  combine  to  produce  a rather  long 
interval  just  before  the  passage  of  the  peak  concentration  when  the  local 
concentration  will  remain  nearly  constant.  In  later  tests  some  local 
concentrations  reached  their  peaks  even  before  the  peak  concentration 
of  the  contaminant  distribution  had  passed. 

The  data  are  insufficient  to  show  conclusively  the  form  of  the 
local  decrease  in  concentration  with  time  for  the  later  tidal  cycles; 
although  other  evidence  indicates  that  the  local  decay  curve  is  not 
exponential  for  large  times. 

Test  4 was  run  under  the  same  conditions  as  Tests  1-A  and  1, 
except  for  river  flow,  which  was  set  at  the  relatively  low  value  of  3000 
sec  ft  at  Trenton.  The  general  features  of  the  contaminant  distribution 
curves,  as  shown  in  Figure  8,  were  similar  to  those  found  for  Tests  1-A 
and  1.  However,  not  only  was  the  rate  of  downstream  movement  decreased 
due  to  the  decrease  in  river  flow,  but  the  rate  of  decrease  in  peak  concen- 
tration is  appreciably  less  in  Test  4 than  for  Test  1=A  or  Test  1.  Figure  9 
shows  the  change  in  peak  concentration  with  time  for  Test  4.  The  rapid 
decrease  in  concentration  during  the  first  tidal  cycle  is  simiia-  that 
which  occurred  in  the  earlier  tests.  The  plot  cf  the  logarithm  of  concen- 
tration versus  the  time  in  tidal  cycles  after  the  second  tidal  cycle  also 
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appears  as  a straight  line,  as  in  the  case  of  Test  1-A  and  Test  1.  However, 
the  slope  of  this  line  in  Test  4 is  only  one-half  the  slope  of  the  decay  curves 
for  the  earlier  tests. 

The  movement  of  the  peak  concentration  with  time  for  Test  4 is 
shown  in  Figure  10.  Again  we  see  that  the  peak  concentration  moves  at 
a slower  rate  than  the  calculated  net  water  movement,  though  the  differences 
are  less  here  than  in  Tests  1-A  and  1. 

The  local  variation  of  contaminant  concentration  with  time  is 
shown  in  Figures  11  and  12  for  selected  stations  above  and  below  the 
release  point.  These  figures  also  reveal  the  slower  rate  of  decrease  in 
the  contaminant  concentration  for  this  test  as  compared  to  Tests  1-A  and  1. 

The  rate  of  river  inflow  thus  appears  to  exert  a significant  in- 
fluence upon  the  rate  of  decrease  in  the  contaminant  concentration. 

Tests  2^  and  3^  - Release  Point  292,  5 

Station  292.5  is  well  into  the  estuary  with  salinities  of  about  15%o 
for  mean  river  flow.  Test  2 was  made  with  mean  tide  and  mean  river 
flow  (12,350  sec  ft  at  Trenton)  while  Test  3 was  made  with  mean  tide  and 
high  river  flow  (31,300  sec  ft  at  Trenton). 

Unfortunately,  much  of  the  data  from  Test  2 is  unusable  since  large 
amounts  of  dye  were  adsorbed  on  the  sample  bottles.  In  addition,  the 
density  of  the  dye  being  less  than  the  density  of  the  water  at  the  release 
point  the  contaminant  was  seen  to  "mushroom"  and  spread  out  in  the  surface 
layer  during  the  earlier  tidal  periods.  (In  Test  3 the  density  of  the  dye 
was  adjusted  to  that  of  the  surrounding  water  and  this  difficulty  was  not 
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repeated.)  Only  the  data  for  tidal  cycles  6,  24,  36,  and  42  were  usable. 

The  longitudinal  distributions  of  the  contaminant  concentration 
for  various  tidal  cycles  for  Test  3 are  shown  in  Figure  13.  The  distribu- 
tions recoverable  from  the  data  for  Test  2 were  similar  to  those  for  Test  3. 

Figure  14  shows  the  movement  of  the  peak  concentration  for  Test  2 
as  well  as  the  computed  mean  downstream  water  movement.  In  Test  2 
where  the  dye  "mushroomed"  into  the  surface  layers,  the  downstream 
movement  of  the  peak  concentration  was  greater  than  or  equal  to  the  com- 
puted net  water  displacement  during  the  first  ten  tidal  cycles.  This  would 
be  expected  since  the  surface  layers  have  an  average  downstream  motion 
greater  than  the  average  for  the  entire  body  of  water.  As  vertical  mixing 
progressed,  the  motion  of  the  peak  concentration  decreased. 

Figure  15  shows  the  movement  of  the  peak  concentration  and  com- 
puted net  water  displacement  for  Test  3.  Here  again  the  movement  of 
the  peak  concentration  is  slower  than  the  calculated  net  water  flow.  The 
difference  for  this  test  was  the  greatest  observed. 

Figure  16  plots  peak  concentration  against  time  for  Test  3.  It 
also  shows  the  usable  data  from  Test  2.  Although  the  data  from  Test  2 
show  more  scatter  than  the  data  from  Test  3,  the  rate  of  decrease  in 
peak  concentration  seems  to  be  about  the  same  for  both  tests.  As  in 
the  previous  tests  the  relation  seems  to  occur  in  two  parts,  an  initial 
rapid  decrease  followed  by  a Blower  exponential  decay. 

Since  so  much  data  were  lost  in  Test  2 a conclusive  evaluation  of 
the  effects  of  changed  river  flow  on  concentration  is  not  possible.  How- 
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ever,  it  seems  that  the  effect  at  Station  292.  5 may  not  be  so  great  as  at 
Station  52.5. 

Curves  showing  the  local  change  in  contaminant  with  time  in  Test  3 
at  stations  above  and  below  the  release  point  are  given  in  Figures  17  and 
18  respectively. 

Tests  f> , 6 , and  7_  ~ Release  Point  111 

These  three  tests  were  conducted  to  determine  the  effects  of  change 
in  tidal  amplitude  on  the  time  change  in  contaminant  concentration.  The 
dye  volume,  with  a concentration  of  1000  ppm,  was  introduced  at  high 
water  slack  at  Station  111  which  is  located  near  Chester,  Pennsylvania. 
This  station  is  very  nearly  the  upper  limit  of  salt  water  intrusion  for  the 
conditions  of  mean  river  flow  under  which  these  tests  were  run. 

Test  5 was  conducted  under  conditions  of  mean  tide,  TeEt  6 under 
conditions  of  spring  tide,  and  Test  7 under  conditions  of  neap  tide.  Since 
the  tidal  velocities  and  tidal  excursions  would  vary  under  these  varying 
tidal  egimes,  it  would  be  expected  that  the  diffusion,  and  hence  the  time 
change  in  concentration,  would  vary  among  these  three  tests.  As  will 
be  pointed  out  below  the  expected  variation  was  neither  clearly  established 
nor  definitely  disproved  by  the  data. 

The  longitudinal  distributions  of  the  contaminant  at  various  tidal 
cycles  for  these  three  tests  are  given  in  Figures  19,  20,  and  21.  These 
distribution  curves  are  similar  to  those  encountered  in  previous  tests  in 
that  they  have  the  general  shape  of  a normal  function,  and  the  change  in 
peak  concentration  with  time  is  large  at  first,  becoming  smaller  for 
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later  tidal  cycles. 

The  downstream  movement  of  the  peak  con  entration  was  essen- 
tially the  same  for  all  three  tests,  and  is  shown  in  Figure  22.  As  with 
all  the  other  tests,  the  observed  movement  of  the  peak  is  at  a slower  rate 
than  the  computed  downstream  water  movement. 

Figures  23,  24,  and  25  are  semilog  plots  of  the  peak  contaminant 
concentration  versus  time  for  each  of  the  three  Tests  5,  6,  and  7 respec- 
tively. Subsequent  to  the  first  two  tidal  cycles,  the  semilog  plots  of  peak 
concentration  versus  time  could  be  adequately  described  by  a straight 
line  for  the  first  five  tests.  For  Tests  5,  6,  and  7,  however,  it  appears 
that  this  simple  exponential  time  variation  is  inadequate.  In  Figure  24 
the  curve  from  Test  1-A  is  also  entered  to  show  the  sharp  difference 
between  the  characteristic  time  change  in  concentration  of  the  first  five 
tests  as  compared  to  these  last  three  tests. 

The  local  time  variation  in  concentration  at  selected  stations  is 
shown  for  these  three  tests  in  Figures  26,  27,  28,  29,  30,  and  31. 

The  plote  of  peak  concentration  versus  time  do  not  show  a systematic 
difference  from  test  to  test  which  can  be  readily  ascribed  to  the  differ- 
ences in  tidal  regimes.  One  hypothesis  concerning  the  effect  of  tide  is 
that  with  higher  tidal  velocities  mixing  would  be  increased  and  hence 
the  peak  concentration  should  decrease  at  a faster  rate,  at  least  during 
the  early  stages.  It  is  difficult  to  measure  accurately  the  actual  peak 
concentration  at  the  first  tidal  cycle,  and  also  after  about  30  tidal  cycles. 
Some  information  is  obtained  by  comparing  the  peak  concentrations  for 


the  three  tests  for  tidal  cycles  after  the  1st  and  prior  to  the  30th.  In 
Table  I the  peak  concentrations  for  the  3rd,  6th,  12th,  and  20th  tidal 
cycles  subsequent  to  release  are  given  for  each  of  the  three  tests. 


Table  I 


Effects  of  Tidal  Regime 


Tidal  Cycle 

Concentration  in  ppm  for  Tests  5, 

6,  and  7 

Test  5 

HP *.  / 

ACDt  U 

Test  7 

Mean  Tide 

Spring  Tide 

Neap  Tide 

3 

5.45 

4.  58 

5.81 

6 

3.25 

2.  83 

3.34 

12 

2.09 

1.42 

1.75 

20 

1.  05 

0.  76 

0.  84 

The  concentrations  for  all  four  tidal  cycles  are  lowest  for  spring 
tide  conditions  when  mixing  should  be  at  a maximum.  The  concentrations 
for  tidal  cycles  3 and  6 are  highest  for  the  neap  tide  conditions,  as  would 
be  expected,  but  for  tidal  cycles  12  and  20  the  concentrations  are  highest 
for  Test  5,  which  was  run  under  mean  tidal  conditions.  The  tidal  mech- 
anism for  the  model  failed  after  cycle  30  during  Test  5,  and  it  may  be 
that  a slight  malfunctioning  occurred  prior  to  the  final  failure,  resulting 
in  somewhat  erroneous  results.  In  any  case  the  data  presented  in  Table  I, 
while  not  conclusive,  do  suggest  that  the  decrease  in  concentration  is  at 
a more  rapid  rate  the  greater  the  tidal  velocities,  at  least  during  the 
early  stages  of  mixing. 
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Summary  Discussion 

The  composite  information  from  all  eight  tests  leads  to  the  follow- 
ing general  statements  concerning  the  distribution  and  rate  of  change  of 
distribution  of  the  dye  contaminant. 

1.  Within  the  first  tidal  cycle  subsequent  to  release,  the  con- 
taminant, which  was  in  all  cases  introduced  at  high  water  slack  as  a 
cylindrical  volume  having  a prototype  diameter  of  740  feet,  spread  over 
a longitudinal  distance  approximately  equal  to  a tidal  excursion  ( about 
40,000  feet)  and  laterally  from  shore  to  shore. 

2.  After  the  first  tidal  cycle  the  longitudinal  distribution  of 
the  contaminant  approximated  a normal  curve. 

3.  Lateral  and  vertical  variations  in  the  contaminant  concentration 
were  observed  during  the  early  stages  of  each  test.  The  sampling  pro- 
cedure did  not  allow  a detailed  investigation  of  the  lateral  or  vertical 
variations  in  concentration,  though  the  pertinent  features  that  were 
noted  are  listed  below: 

a.  During  the  early  stages  of  mixing,  relatively  concen- 
trated segments  of  the  contaminant  volume  were  observed  to  become 
entrapped  by  eddies  within  small  embayments  or  behind  projecting  marine 
structures.  These  entrapped  contaminant  segments  would  then  feed  out 
contaminant  into  the  main  stream  after  the  main  mixing  volume  had  been 
carried  out  of  the  immediate  area  by  the  tidal  currents.  It  appeared 
that  this  feature  played  a major  role  in  contributing  to  the  horizontal 


spread  of  the  contaminant. 
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b.  At  prominent  bends  in  the  river  the  advanced  portion  of 
the  contaminant  volume  was  observed  to  cling  closer  to  that  shore  having 
a more  rapid  downstream  movement  of  water. 

c.  In  Test  3,  which  was  conducted  in  that  part  of  the  estuary 
which  exhibits  a net  two-layered  flow  pattern  similar  to  that  described 

by  Pritchard  {1952]'  for  the  Chesapeake  Bay,  the  upstream  spread  of  the 
contaminant  was  observed  to  occur  primarily  along  the  bottom.  During 
the  early  stages  of  the  tests  the  vertical  distribution  for  stations  above 
the  release  point  ehowed  increasing  concentration  with  depth,  while  those 
below  the  release  point  in  general  showed  decreasing  concentration  with 
depth.  After  a few  tidal  cycles,  vertical  mixing  had  destroyed  these 
general  vertical  variations  in  contaminant  concentration. 

4.  The  change  in  peak  concentration  is  characterized  by  a very 
rapid  rate  of  decrease  during  the  first  tidal  cycle  followed  by  a less 

rapid  rate  of  decrease.  In  the  first  five  tests,  namely  1-A,  1,  2,  3,  and  4, 
the  variation  in  peak  concentration  subsequent  to  the  first  tidal  cycle 
can  be  represented  by  a simple  function  of  time.  In  Tests  5,  6,  and  7, 
however,  the  variation  of  peak  concentration  subsequent  to  the  first  tidal 
cycle  cannot  be  adequately  represented  by  a simple  exponential  rate.  The 
reason  for  the  difference  between  these  groups  of  tests  is  not  at  present 
evident. 

5.  The  downstream  movement  of  the  peak  concentration  is  sig- 
nificantly slower  than  the  calculated  velocity  of  the  mean  downstream 


we  ter  movement. 
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Theoretical  Considerations 

A Preliminary  Analysis  of  the  Characteristic  Features  of  the  Flushing 
Study  in  the  Delaware  Model 

The  semilog  plots  of  peak  concentration  versus  time  in.  tidal  cycles 
for  Tests  1-A,  1,  2,  3,  and  4 (Figures  4,  9,  and  16)  are  all  of  the  form 
shown  in  the  following  diagram: 


\ 


Since  this  is  a semilog  plot  the  straight  lines  indicated  by  (A)  and  (B) 
maybe  expressed  formally  using  the  slope-intercept  form  of  the  straight 
line  as 


and 


108 

108  Xh 


= “b|t  + log  3ii 
= -bft  + log  af 


re  spectively 


where  -b^  , 


-b£  are  slopes  and  log  a^,  log  a^  are  the  intercepts.  Since 
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D 

log  A = B implies  ec  = A,  we  may  write: 


\ -bit 

a - aie 

Hnd  = afe=t>ft 


If  we  sum  these  two  equations 


for  the  initial  period 
for  the  final  period 


i.  + 


-b-t 
= a^e  i 


a£e 


•bft 


we  may  use  the  result  as  an  approximation  to  as  indicated  by  the 

full  line  since  to  the  left  of  the  transition  zone  (C)  the  contribution  of 
(B)  to  (A)  is  negligible  and  conversely.  Within  the  transition  zone,  (C), 
(A),  and  (B)  will  both  contribute  and  the  sharp  corner  will  be  rounded  off. 

Th*»T*»fore , formally  at  least,  peak  concentration  as  a function  of 
time  may  be  written 


(i) 


^ = aje-V 


, -bft 

+ a£e  1 


As  pointed  out  above,  a peculiarity  of  this  function  is  that  for 
values  of_t_  outside  of  the  neighborhood  of  t = 1 only  one  of  the  terms 
makes  a major  contribution  to  the  value  of  . Thus  during  the  first 

tidal  cycle  subsequent  to  release,  the  very  rapid  rate  of  decrease  charac- 
terized by  the  first  term  of  equation  (1)  predominates,  and  we  have 


(2) 


= a,e~V 


(first  tidal  cycle  only) 


After  the  second  tidal  cycle,  this  first  term  becomes  unimportant, 
and  the  change  in  peak  concentration  is  characterized  by  the  final  rate 


of  decrease  given  by 
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(3) 


= a. 


rhft 


(after  the  second  tidal  cycle) 


The  variations  in  peak  concentration  for  Tests  5,  6,  and  7,  sub- 
sequent to  the  second  tidal  cycle,  do  not  appear  to  be  adequately  des- 
cribed by  the  exponential  expression  (3).  For  all  tests,  however,  it 
appears  that  the  expression 


(4) 


r 


•* 


= b' 


holds.  Here  b1  and  n are  constants.  The  form  of  the  solution  of  this 
differential  equation  depends  on  the  value  of  n.  For  the  special  case 
of  n = -1,  the  solution  of  {4)  is 


(5) 


in  t = a'  + b't  or  X = 


ae 


b't 


which,  when  b:  is  identified  with  -la  , is  seen  to  be  the  expression  which 
satisfied  the  observed  concentration-time  variations  for  Tests  1-A,  1,  Z, 
3,  and  4. 

For  all  n other  than  n = -1,  the  formal  solution  of  (4)  is 

b+l 


(6) 


= bt  + a 


b = b'  (n+ 1) 


Designating  (n  + 1)  by  -m,  this  expre  ssion  becomes 


(?) 


X 


■ m 


= bt  + a 


By  solving  for_t,  (7)  can  be  expressed  as 


t = (i/b)  y 


V -m 


f 


+ C 


(8) 


C = -a/b 
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or 


In  (t  - c)  = 


"A 


-lnb 


Thus  a plot  of  ln^J?^  versus  In  (t-c)  should  produce  a straight 
line,  of  slope  m.  It  is  possible  to  obtain  an  approximate  value  of  the 
constant  £ graphically  from  a plot  of  data  which  is  suspected  of  showing 
a functional  variation  given  by  (8).  This  has  been  done  for  Tests  5,  6, 
and  7,  and  the  resulting  log-log  plots  of  the  concentration  versus  the  sum 
of  the  time  plus  the  constant  £ are  given  in  Figures  32,  33,  and  34.  It  is 
seen  that  the  relationship  is  reasonably  linear,  and  the  equations  for 
the  three  sets  of  data  are: 


Test  5 
Test  6 
Test  7 


■ 0.637 


--0.914 


•0. 740 


= 0.  0341  t + 0.  225 
= 0.0513  t + 0.0765 
= 0.0506  t + 0.  114 


Since  in  equation  (7)  n = -m  °1  , n varies  between  -1.63  and  -1.92 
for  Tests  5,  6,  and  7,  while  for  the  five  previous  tests  we  found  n = -1. 

For  the  present  we  will  confine  ourselves  to  a discussion  of  the 
characteristic  features  of  the  decay  expression 


(9a) 


which  can  be  expressed  as 


(9b) 


In 


bt 


a'  = In  a 


A plot  of  In  against  time  would  result  in  a straight  line  of 

slope  -b  . The  decay  of  radioactive  material  is  represented  by  a function 
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of  this  type,  and  borrowing  from  the  nomenclature  of  the  nuclear  field  we 
may  consider  a characteristic  time  related  to  the  rate  of  flushing  of  a con- 
taminant which  we  will  call  the  "flushing  half-life.  " The  flushing  half-life 
is  defined  as  the  time  interval  required  to  reduce  a given  concentration 
to  a concentration  one-half  as  great.  It  should  be  noticed  that  if  the  rate 
of  change  of  concentration  is  given  by  equation  (9),  then  the  half-life  is 
a characteristic  parameter  independent  of  time  or  concentration;  that  is, 
if  ^5^  represents  the  concentration  at  time  tp  the  same  interval  of 
time  *s  re<luiretl  to  reduce  the  concentration  to  l/2  ( )>  as  is 


reauired  to  reduce  the  concentration  from  l/2  ( ),  to  l/4  ( ), 

^T»  ' T' 

from  1/4  ( ),  to  1/8  ( ^ ). 

The  half-life  is  readily  obtained  from  the  slope  b,  since 


or 


(10) 


^ tl/2  In  (1/2) 


The  flushing  half-life  is  a convenient  parameter  for  discussing 
certain  features  of  the  flushing  tests.  In  Table  II  the  half-life  for  the 
various  tests  and  for  the  various  characteristic  nortions  of  each  test  is 
given.  It  is  seen  that  the  half-life  characteristic  of  the  rate  of  concen- 
tration decrease  during  the  first  tidal  cycle  is  only  of  the  order  of  l/lOO 
of  the  half-life  associated  with  the  rate  of  concentration  decrease  sub- 
sequent to  the  second  tidal  cycle.  Also  we  note  that  the  half-life  for 
Test  4,  which  was  run  under  low  river  conditions,  is  two  times  the 
corresponding  half-life  for  Tests  1-A  and  1,  which  were  run  under  con- 


- 19 


ditions  of  mean  river  flow,  other  conditions  being  the  same  for  the  three 
tests. 


Table  II 


Test 

Half-life  in 

tidal  cycles 

For  first 
tidal  cycle 

After  second 
tidal  cycle 

1-A  and  1 

0.  16 

10.6 

2 and  3 

0.  11 

13.5 

4 

0.  15 

24.4 

5* 

0.  16 

11.  1 

6* 

0.  15 

10.  8 

7* 

0.  16 

9.8 

30C 

Tests  5,  6,  and  7 are  included  here  for  comparison,  even 
though  the  peak  concentration  could  not  be  adequately  described 
by  a simple  exponential  function  of  time.  The  calculations 
given  are  based  on  the  concentration  change  from  the  fourth 
to  the  fortieth  tidal  cycles,  neglecting  the  departure  from  the 
exponential  form  between  these  points. 


Some  Theoretical  C onsiderations  Suggested  by  the  Model  Study  of  Flushing 
It  can  be  shown  (Pritchard,  1952,  unpublished  notes)  that  the  one- 
dimensional expression  for  the  time  rate  of  change  of  mean  concentration 
in  a channel  of  varying  cross  section  may  be  expressed  as 

(13)  ' 'vlf  + ]\  z*lAK/ ■ Ixl 

where  x is  the  coordinate  axis  along  the  axis  of  the  channel,  v the  mean 
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velocity  in  the  x-direction,  A the  cross  sectional  area,  and  K, 


the 


horizontal  eddy  diffusivity.  If  it  is  assumed  that  the  contaminant  dis~ 
tribution  approximates  the  normal  curve,  then  it  can  also  be  shown  that 


(14) 


K 


- O' 


where  is  the  peak  concentration  of  the  contaminant. 

In  Tests  1-A,  1,  2,  3,  and  4 it  was  found  that 


In 


a**bt 


equation  (4) 


Therefore,  for  these  tests. 


(15) 


K 


b c r 


and  the  eddy  diffusivity  varies  as  the  square  of  the  standard  spread 

of  the  contaminant  distribution. 

In  Tests  5,  6,  and  7,  however,  it  was  found  that 


X 


."H 


a + bt 


equation  (7) 


for  which  the  diffusivity  would  be  given  by 


(16) 


Kr  ~ where  /£ 


b / M \m 

■sr  ( ~J^r  ’ 


A^  is  the  mean  cross  sectional  area  of  the  segment  of  the  estuary  occupied 
by  the  contaminant  volume,  and  M_  is  the  total  weight  of  the  contaminant. 
From  equation  (16)  it  can  be  seen  that,  for  these  tests,  the  diffusivity 
varies  inversely  as  the  mth  power  of  the  cross  sectional  area  and 
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directly  as  the  (2-m)th  power  of  the  standard  spread  of  the  contaminant 
distribution. 

In  Tests  5,  6,  and  7 m had  an  approximate  range  from  0.5  to 
1.0  showing  extremes  for  the  diffusivity  of 


(17a) 

(17b) 


■#s<r 
* %* 


i.g 


and 


Stommel  (1949)  has  suggested  that  if  L represents  a scale  factor 
characterizing  the  size  of  the  area  within  which  eddy  diffusion  is  being 
investigated,  then  the  eddy  diffusivity  is  approximately  proportional  to 
the  4/3  power  of  L_.  This  concept  was  developed  for  diffusion  in  an  area 
unrestricted  by  lateral  boundaries  and  would  hence  seem  inapplicable  to 
this  study.  It  is  interesting-  therefore,  to  note  that  the  mean  exponent  of 
for  Tests  5,  6,  and  7 is  approximately  equal  to  the  4/3  power 
suggested  by  Stommel. 

That  the  movement  of  peak  concentration  was  in  all  cases  slower 
than  the  calculated  net  water  movement  calls  for  comment.  Pritchard 
(1952,  unpublished  notes)  has  shown  that,  if  the  distribution  of  the  con- 
taminant approximates  a normal  curve,  then  the  position  of  the  peak 
concentration  is  given  approximately  by 

^ M1  S' ! J ^ftu  1 

U8,  YW  = ' £^“35" 


m is  the  position  of  the  peak  concentration  when  the  origin  x = 0 is  set 
at  the  position  of  the  peak  concentration  at  time  t = 0.  is  the  net 


# 
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movement  due  to  the  velocity  of  the  water,  i.e.  D = J vdt  where  v 

‘'o 


rl 

, i.e.  D = I vdi 
'■'o 

is  the  mean  water  velocity.  The  subscript  "o"  indicates  the  value  of 
the  parameter  at  time  zero.  The  meanings  of  the  other  symbols  are  as 
previously  explained. 

In  this  study  the  factor  in  braces  in  the  second  term  of  the  right- 
hand  number  of  equation  (18)  was  always  positive  and  therefore  rr>  <D* 
is  to  be  expected. 

The  development  of  a suitable  theoretical  description  of  other 
aspects  of  these  flushing  tests  is  in  progress.  However,  in  view  of  the 
unexplained  difference  between  the  characteristic  rates  of  decay  for  the 
first  five  tests  and  for  the  last  three  tests,  further  experiments  with  the 
model  would  be  desirable. 


Suggestions  Regarding  Flushing  in  the  Prototype 

The  question  now  arises  as  to  the  extent  to  which  the  results  ob- 
tained from  the  model  might  be  applied  to  the  prototype.  In  this  regard 
we  may  consider  certain  features  of  the  model  verification. 

The  model  was  constructed  according  to  Froude  scaling.  The  model 
velocity  structure  and  salt  distribution  which  would  result  from  Froude 
scaling  alone  would  not  be  comparable  to  the  prototype,  since  the  Reynolds 
scale  would  be  far  from  satisfied.  Tn  order  to  increase  the  turbulence, 
the  model  was  artificially  roughened.  The  extent  of  this  artificial  rough- 
ening was  empirically  determined  by  roughening  until  the  velocity  distri- 
bution and  the  tidal  heights  matched  those  known  to  exist  in  the  prototype. 


- 23 


It  was  then  found  that  the  salinity  distribution  in  the  model  satisfactorily 
matched  that  found  in  the  prototype  under  various  river  flow  regimes. 

This  salinity  verification  implies  that  the  artificial  roughening  of 
the  model  successfully  sealed  those  eddies  in  the  turbulent  regime  which 
are  important  in  maintaining  the  salt  distributions.  There  is  strong  indi- 
cation, however,  that  the  scale  of  the  turbulent  eddy  system  which  is 
involved  in  any  particular  mixing  process  depends  on  the  size  of  the  con- 
taminant volume  being  mixed.  Therefore,  the  eddies  which  would  be 
involved  in  the  mixing  of  a highly  concentrated  contaminant  spread  over 
only  a small  area  may  not  be  properly  scaled  in  the  model,  even  though 
the  eddies  important  in  maintaining  the  salt  distribution  are  properly 
scaled  from  the  prototype  to  the  model. 

On  the  basis  of  the  above  reasoning,  the  following  preliminary  con- 
clusions regarding  the  application  of  the  results  of  the  model  studies  on 
flushing  to  prototype  conditions  are  presented. 

1.  It  is  doubtful  that  the  extremely  rapid  rate  of  concentration 
decreases  observed  during  the  first  tidal  cycle  in  the  model  studies 
actually  represent  in  magnitude  the  corresponding  initial  rate  of  mixing 
which  would  occur  in  the  prototype.  Immediately  after  introduction  of 
the  dye,  large  gradients  exist  and  molecular  diffusion  would  be  of  some 
importance  in  the  model,  though  less  so  in  the  prototype  since  there  can 
be  no  scaling  of  this  phenomenon  from  the  prototype  to  the  model.  Like- 
wise, as  discussed  above,  the  eddy  scale  important  to  the  mixing  of  the 
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restricted  contaminant  volume  during  this  early  phase  may  not  be  properly 
scaled  in  the  model. 

Undoubtedly,  the  rate  of  concentration  decrease  during  the  early 
phase  of  mixing  in  the  prototype  will  be  substantially  greater  than  at 
later  stages,  since  the  initial  contaminant  volume  may  spread  laterally 
as  well  as  longitudinally  (which  is  the  case  in  the  model  also).  However, 
the  rapid  decrease  in  concentration  in  the  model  during  the  first  tidal 
cycle  is  substantially  greater  than  would  be  expected  as  a result  of  the 
coupling  of  the  lateral  spread  with  the  longitudinal  mixing.  This  is  a 
feature  which  may  allow  prototype  verification,  since  it  is  possible  that 
tracer  material  which  could  be  followed  for  one  or  two  tidal  cycles  might 
be  safely  added  to  the  natural  estuary. 

2.  The  rate  of  concentration  change  observed  in  the  model  sub- 
sequent to  the  second  tidal  cycle  is  probably  applicable  to  the  prototype. 

By  this  time  the  contaminant  had  spread  from  shore  to  shore  laterally 
and  over  approximately  a tidal  excursion  longitudinally.  The  eddy  regime 
contributing  to  the  mixing  would  be  of  about  the  same  order  of  magnitude 
as  that  which  is  important  to  the  maintenance  of  the  salt  distribution,  and 
hence  probably  properly  scaled  from  prototype  to  model. 

The  flushing  half-lives  given  in  Table  II  for  the  period  subsequent 
to  the  second  tidal  cycle  after  dye  introduction  are  therefore  probably 
representative  of  the  expected  rates  of  flushing  in  the  prototype  for  con- 
taminant distributions  which  are  spread  over  approximately  a tidal  excursion. 

3.  The  suggestion  that  the  eddy  ciffusivity  varies  with  at  least  the 


« 
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first  power,  and  possibly  the  second  power,  of  the  standard  spread  of  the 
contaminant  distribution  throws  some  doubt  on  the  application  of  techniques 
for  computing  flushing  characteristics  of  an  estuary  which  are  based  on 
the  salinity  distribution,  or  conversely,  on  the  observed  or  computed 
distribution  of  freshwater.  Though  these  techniques  may  adequately 
describe  the  fresh  water-salt  water  distribution,  they  may  not  describe 
the  concentration  change  of  a contaminant  which  is  introduced  into  the 
estuary  with  an  initial  distribution  considerably  different  from  the  existing 
fresh  water  or  salt  water  distribution. 
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Figure  1 


Location  Map  for  Model  Study  of  the  Delaware  Estuary. 


) 

Figure  2.  Schematic  representation  of  dye  entrappment  by  shore  features. 

Concentration  is  indicated  by  intensity  of  shading.  As  peak  concentration  of 
main  dye  volume  moves  downstream  on  the  ebb  tide,  a portion  of  the  dye 
spreads  into  shore  indenture  _A,  as  in  diagram  (i).  This  portion  of  the  dye 
is  entrapped  by  the  indenture  as  the  main  dye  volume  moves  on  downstream, 
as  shown  in  diagram  {ii).  In  diagram  (lii)  the  entrapped  dye  feeds  out  into 
the  main  channel,  this  contributing  to  the  longitudinal  spread.  The  main  dye 
volume  is  shown  in  this  diagram  to  have  spread  into  a second  shore  indenture. 

In  diagram  (iv)  the  phenomena  is  shown  repeated  for  a second  indenture  B on 
the  flooding  tide. 
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APPENDIX  I 


Description  of  the  Dye  Release  Gear,  the  Water  Sampling  Equipment, 
the  Method  of  Sampling,  and  the  Determination  of  Dye  Concentration 

Dye  Release  Gear 

The  gear  used  to  :,elease  the  methylene  blue  used  in  these  experi- 
ments is  shown  in  Figure  A-I-l.  It  consists  of  an  open  brass  cylinder 
held  firmly  against  a rubber  lined  base  by  means  of  a two  inch  wide  bar 
and  two  outrigged  rods.  The  rods  were  screwed  into  threaded  holes  in 
projections  of  the  bottom  piece,  and  pulled  down  on  the  two  inch  bar 
which  was  placed  across  the  top  of  the  cylinder.  The  bottom  of  the 
cylinder  was  thus  drawn,  tightly  down  on.  the  rubber  sheet  which  was 
bonded  to  the  base  and  leakage  from  the  container  was  prevented. 

At  the  appropriate  release  point  a recess  was  made  in  the  bottom 
of  the  model  of  just  the  correct  size  and  shape  to  hold  the  bottom  piece 
of  the  container.  The  depth  of  the  water  at  the  release  point  was  accurately 
measured  at  high  water,  and  the  height  of  dye  in  the  container  was  made 
to  agree  with  this  depth  of  water. 

During  the  flooding  tide  just  prior  to  the  release  time  the  dye 
container  was  introduced  into  the  model,  the  bottom  piece  resting  in 
the  recession.  The  rods  were  unscrewed  from  the  bottom  and  withdrawn, 
while  one  man  held  his  foot  securely  on  the  rod  across  the  top  of  the 
container,  bearing  down  with  sufficient  force  to  prevent  leakage  of  dye 
from  the  bottom  of  the  container.  At  precisely  high  water  slack  cf  the 
tidal  period  chosen  for  release,  the  cylinder,  now  no  longer  secured  to 
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the  bottom  piece,  was  lifted  quickly  but  carefully  upward  from  the  model, 
leaving  the  bottom  piece  remaining  in  the  recess  in  the  model  bottom, 
and  the  cylindrical  dye  volume  suspended  in  the  water  of  the  model.  Since 
the  recess  had  been  made  to  just  fit  the  bottom  piece  of  the  container, 
the  bottom  of  the  model  with  the  bottom  piece  in  place  was  without  any 
indentation  or  protrusion  resulting  from  the  release  of  the  dye. 

Sampling  Equipment 

Since  the  model  has  a vertical  scale  of  1 to  100,  0.0.1  foot  of  depth 
in  the  model  corresponds  to  1 foot  of  depth  in  the  prototype.  The  sampling 
equipment,  therefore,  had  to  be  made  relatively  small  so  as  not  to  inter- 
fere with  the  flow. 

The  basic  unit  of  the  sampling  gear  was  a thin  glass  tube,  bent  at  a 
right  angle  near  the  end  and  drawn  to  an  elongated  point.  The  glass  tubing 
was  connected  by  plastic  tubing  to  a 50  ml  syringe.  The  assembly  could 
be  operated  by  hand  as  a single  unit,  or  grouped  together  with  other  simi- 
lar assemblies  to  give  a multiple  depth  sampler.  In  any  case  the  glass 
tubing  was  taped  to  a metal  rod  so  that  the  metal  rod  projected  at  least  a 
slight  amount  below  the  lowest  glass  tip. 

In  the  single  unit  hand  samplers,  the  glass  tube  was  taped  to  the 
metal  rod  so  that  the  bent  tip  was  about  0.  2 ft  from  the  bottom  of  the 
rod  corresponding  to  a prototype  distance  above  tr. e bottom  of  20  ft.  In 
this  manner  the  rod  could  be  placed  with  its  end  on  the  bottom  at  the 
desired  sampling  position  ir.  mid  channel  ar.d  the  sample  would  be  drawn 
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from  about  middepth,  since  the  channel  up  to  Philadelphia  is  maintained 
at  the  prototype  depth  of  40  feet. 

The  multi-depth  sampling  gear  is  shown  in  Figure  A- 1-2.  The 
glass  tubes  were  taped  to  a metal  rod  so  that  the  bent  tips  occurred  at 
appropriate  intervals  --  usually  surface,  middepth,  and  bottom.  The 
rod  was  part  of  a tripod  assembly  usually  employed  in  measuring  the  sur- 
face elevation  of  the  water  in  the  model,  and  the  sampling  assembly  was 
thus  supported  in  position  by  the  tripod.  The  plastic  tubes  leading  from 
the  three  glass  tubes  were  connected  to  three  50  ml  syringes  secured  to  a 
wooden  rack.  The  pistons  for  the  syringes  were  all  clamped  to  a wooden 
bar,  which  in  turn  was  connected  to  a cable  leading  over  a pulley  at  the 
top  of  the  wooden  rack  . This  cable  was  attached  to  a bar  secured  to 
three  household  rattraps. 

A spring  latch  held  the  three  rattraps,  which  were  connected  in  the 
open  or  "set'1  position.  In  this  position  the  cable  was  just  taut  with  the 
pistons  in  the  syringe  all  the  way  in.  When  the  latch  was  released,  the 
spring  on  the  rattraps  would  pull  down  on  the  cable,  and,  hence,  up  on 
the  syringes  and  the  sample  would  be  drawn. 

In  Test  1-A,  use  was  made  of  a motor  driven  multi~sampler  which 
took  eleven  simultaneous  samples,  drawing  the  sample  evenly  over  the 
7-1/2  minute  tidal  period  of  the  model.  This  mechanism  was  not  employed 
after  the  first  test,  since  there  was  some  difficulty  in  fitting  the  results 
from  these  samples  with  the  samples  taken  at  high  water  slack. 


- 64  - 


Sampling  Procedure 

Before  introduction  of  the  dye,  the  rattrap  samplers  were  set  up 
at  appropriate  stations  above  and  below  the  release  point.  Other  locations 
were  chosen  where  single  hand-drawn  samples  were  to  be  obtained.  On 
selected  high  water  slack  periods,  subsequent  to  the  release,  samples 
were  collected  at  the  designated  positions. 

Since  the  tide  in  the  Delaware  is  partially  progressive,  a particular 
phase  of  the  tide  occurs  downstream  first  and  can  be  followed  up  the  Bay  and 
the  River.  Small  strips  of  paper  were  strewn  on  the  surface  of  the  model, 
and  the  time  of  high  water  slack  determined  at  each  particular  station  by 
noting  he  time  when  these  floats  stopped  moving  upstream.  Thus,  the 
sampling  was  not  simultaneous  at  all  stations,  but  related  to  the  tidal 
stage,  the  first  samples  being  taken  at  the  downstream  edge  of  the  spread 
of  the  dye,  and  the  last  samples  for  any  particular  tidal  cycle  at  the 
upstream  end  of  the  spread.  With  several  individuals  participating  in 
the  sampling,  it  was  quite  possible  to  keep  pace  with  the  progressing 
tidal  phase  as  it  moved  upstream. 

The  samples  were  transferred  from  the  syringes  into  wide  mouth 
quart  jars.  As  will  be  pointed  out  in  a later  appendix,  this  was  far 
from  the  most  satisfactory  procedure.  The  jars  were  then  assembled 
for  analysis. 

Determination  of  the  Dye  Concentration 

The  concentration  of  dye  was  determined  photometrically  using  both 
a Beckman  model  DU  spectrophotometer  and  a Coleman  spectrophotometer. 
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At  the  beginning  of  each  run  a sample  was  taken  from  the  dye  volume  to  be 
introduced,  and  appropriate  volumetric  dilutions  were  made  from  this 
sample.  The  Beers-law  calibration  curve  of  adsorption  versus  concen- 
tration was  then  obtained  for  each  instrument,  using  the  diluted  samples. 

Samples  of  water  from  the  model  were  also  drawn  prior  to  the  dye 
introduction  to  serve  as  blanks  in  the  calibration. 

The  ground  around  the  model  is  quite  dusty,  and  even  though  the 
model  is  well  washed  down  between  tests,  there  occurs  in  the  model  water 
some  suspended  material  which  appeared  to  interfere  with  the  determina- 
tion of  low  concentrations.  The  procedure  employed  in  eliminating  this 
difficulty  is  essentially  the  procedure  described  by  Carritt  and  Carpenter 
(1951)  in  their  paper  on  corrections  for  turbidity  interference  in  chloro- 
phyll analysis. 

The  procedure  depends  upon  the  contrast  between  the  extinction 
curve  due  to  the  dye  and  the  extinction  curve  resulting  from  the  turbidity. 
This  difference  is  shown  in  Figure  A-I-3.  The  broken  line  A represents 
schematically  the  extinction  curve  for  the  dye  in  distilled  water.  This 
curve  is  characterized  by  a sharp  maximum  in  the  extinction  at  a wave 
length  of  665  xryJU  . The  dashed  line  _B_  is  representative  of  the 
typical  extinction  curve  resulting  from  turbidity.  This  curve  shows  a 
steady  decrease  with  wave  length,  with  no  characteristic  adsorption  peak. 

The  addition  of  these  two  curves  gives  curve  C , drawn  as  a.  con- 
tinuous line.  The  typical  extinction  curve  for  the  dye  is  moved  up  and 
modified  somewhat  by  the  turbidity. 
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Normally  the  extinction  Em  at  the  characteristic  wave  length  /V 
is  taken  as  proportional  to  the  concentration.  However,  if  the  sample  is 
somewhat  turbid,  the  resulting  extinction  E 1 would  indicate  a larger  dye 
concentration  than  actually  existed.  If  two  points,  E&  and  E^  are  taken 
on  the  dye  extinction  curve  well  away  from  the  peak  extinction  at  wave 
lengths  /V  and  respectively,  and  E is  the  linear  interpolation 


between  E and  E at  wave  length 
j*.  b 


* 


m 


then  the  distance  E to  E is 


m 


also  proportional  to  the  concentration.  However,  the  distance  Em-  E^  is 
approximately  equal  to  the  distance  E'm  = E1  on  the  combined  extinction 
curve  for  ^oth  dye  and  turbidity. 

The  procedure  calls  for  running  a calibration  curve  of  the  difference, 

E - E^  versus  concentration.  To  do  this,  measurements  are  made  of  the 

m c 

extinction  readings  at  the  three  wave  lengths  A , A , and  ^A 

a ' 'm  b 

From  the  extinctions  E and  E,  the  correction  E can  be  obtained  since 

a b c 


Ec  = Ea 


<Ea-Eb> 


^m 


A 


A 
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The  difference  Em  - Ec  is  then  taken  as  proportional  to  the  concentration. 
The  unknown  samples  are  analyzed  in  the  same  manner,  and  hence  the 
influence  of  the  turbidity  is  essentially  eliminated. 


Extinction 
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■wave  length 
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Figure  A-I-3.  Effect  of  turbidity  (curve  B)  on 
the  extinction  of  dye.  The  curve  for  the  dye  in 
distilled  water  is  A_  , and  in  water  of  turbidity 
given  by  curve  13  the  combined  extinction  curve 
is  given  by  C . The  distance  (Em-  E ) is  , how  - 
ever,  approximately  equal  to  (E'm-  E^.'  ) , and 
may  also  be  expressed  as  a function  of  dye 
concentration. 
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APPENDIX  11 


Visual  Observations  o£  the  Spread  of  the  Dye. 

Tables  of  Observed  Dye  Concentrations. 

Fox  Tests  3 through  7 the  location  of  the  visible  limits  of  the  dye 
was  noted  and  charts  were  made  during  the  tests.  Notes  were  also  made 
of  any  interesting  features  that  were  observed.  These  charts  are  presented 
in  this  appendix. 

It  was  found  that  the  limits  of  the  dye  concentration  as  ue  te  i- liiiue  u 
visually  corresponded  to  an  average  concentration  of  less  than  0.  1 ppm. 

This  corresponds  to  more  than  a ten  thousand  fold  dilution  of  the  initial 
concentration. 

This  appendix  also  contains  tables  of  the  observed  dye  concentrations. 
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065 
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I 


Station 

Depth 

Cone . 

Station 

Depth 

— 

Cone . 

.(2Pml 

(ppm) 

100 

M(b) 

.008 

100 

M(b) 

.000 

110 

M(b) 

.Ooo 

110 

M(b) 

.001 

130 

S(b) 

.000 

130 

S 

n.d. 

M(b) 

.055 

M(b) 

. 000 

B(b) 

.001 

B(b) 

.000 

150 

M(b) 

.000 

150 

M 

n.d. 

B(b) 

.000 

160 

S 

n.d. 

160 

s 

n.d. 

M(b) 

.009 

M 

n.d. 

B 

n.d. 

B(b) 

.003 

TIDAL 

CYCLE 

46 

180 

S 

n.d. 

180 

S(b) 
M(b ) 

.000 

.008 

M(b) 

B 

.003 

n.d. 

B(b) 

.004 

200 

M(c) 

.026 

200 

S 

n.d. 

M 

n.d. 

220 

S 

n.d. 

B(c) 

.070 

M 

n.d. 

B 

n.d. 

220 

S(c) 

.110 

.014 

M(c) 

.060 

240 

M(b) 

TIDAL 

CYCLE 

B(c) 

.060 

B(b) 

.015 

.045 

40 

240 

S 

n.d. 

260 

M(b) 

M(c) 

.o4i 

B(c) 

. 000 

280 

M(b) 

.001 

260 

S 

n.d. 

300 

S(b) 

.020 

M(c) 

.001 

M(b) 

.000 

B 

n.d. 

B 

.000 

280 

S 

n.d. 

320 

S 

n.d. 

M 

n.d. 

B 

n.d. 

350 

M(b) 

.001 

B 

n.d. 

300 

S(c) 

.04l 

M 

n.d. 

B(c) 

.013 

320 

S 

n.d. 

M 

n.d. 

B 

n.d. 

350 


M 


Ti  • d . 
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TEST  1 


Station 

■g 

Hi 

100 

M(b) 

.Ooo 

110 

M(b) 

.000 

130 

S(b) 

.000 

M(b) 

.000 

B(b) 

.000 

150 

M(b) 

.000 

160 

S(b) 

.000 

M(b) 

.000 

B(b) 

.000 

180 

S 

n.d. 

M(b) 

.000 

B 

n.d. 

200 

S(b) 

.003 

B 

.003 

TIDAL 

CYC  IE 

220 

S(b) 

.009 

52 

M(b) 

.017 

B(b) 

.042 

24o 

S(b) 

.054 

M(b) 

.019 

260 

S(b) 

.010 

M(b) 

.008 

B(b) 

.000 

280 

S(b) 

.022 

M(b) 

.026 

B(b) 

.016 

300 

S(b) 

.010 

M(b) 

.007 

B(b) 

.014 

320 

S(b) 

.005 

M(b) 

.00  6 

350 

S 

.007 

M 

.000 

B 

.000 

Station 

Depth 

EH 

TIDAL 

300 

M 

.045 

CYCIE 

56 

320 

M 

.031 

200 

M 

itn 

O 

0 

220 

M 

.022 

TIDAL 

240 

M 

.044 

CYCLE 

58 

260 

M 

On 

O 

• 

280 

M 

.035 

350 

M 

.023 

I 
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TEST  3 
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TIDAL 

CYCLE 

20 


Station 

Depth 

Cone . 

(ppm). 

220 

M 

.017 

240 

S 

.074 

M 

.041 

B 

.043 

280 

S 

.185 

M 

.250 

B 

.279 

300 

M 

.484 

320 

S 

.365 

M 

.400 

B 

.4io 

320E 

M 

.471 

320W 

M 

.409 

330 

M 

.510 

350 

s 

.360 

M 

.489 

B 

.435 

350E 

M 

.521 

350W 

M 

• 571 

390 

M 

.233 

390E 

M 

.280 

390W 

M 

.301 

430 

s 

.099 

M 

.068 

B 

.070 

430E 

M 

.046 

430W 

M 

.020 

TEST  4 


- 99  - 


Station 

Depth 

Cone . 

(PPm), 

Station 

Depth 

r 

Cone . 
(ppm) 

-30 

M 

Ml 

-50 

M 

.056 

-10 

M 

-30 

M 

.073 

10 

M 

.311 

-10 

M 

.130 

20 

S 

.251 

10 

M 

.230 

M 

.271 

B 

.193 

20 

M 

.166 

40 

S 

.540 

30 

M 

0 

0 

0 

M 

.468 

B 

.389 

TIDAL 

4o 

n 

O 

.283 

CYCU3 

M 

.259 

TIDAL 

60 

S 

.800 

30 

B 

.292 

CYCLE 

M 

• 730 

20 

B 

.745 

60 

S 

.591 

M 

.485 

80 

3 

1.35 

3 

.440 

M 

1.36 

B 

1.20 

80 

S 

.84o 

M 

.800 

100 

S 

1.36 

B 

.84-0 

M 

1.2c 

B 

1.29 

100 

S 

.943 

M 

.928 

120 

S 

• 790 

B 

1.04 

M 

.705 

B 

.699 

120 

M 

.741 

l4o 

S 

• 331 

l4o 

S 

.588 

M 

.352 

M 

.507 

B 

.337 

B 

.534 

160 

S 

• 315 

160 

M 

.263 

M 

.126 

B 

.126 

180 

S 

.214 

M 

.190 

180 

S 

.085 

B 

.170 

M 

.053 

B 

.081 

200 

M 

.113 

200 

M 

.065 

220 

M 

.101 

240 

M 

.073 

260 


M 


073 

O65 


I 


i 
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TEST  6 


80 

M 

.023 

90 

S 

.000 

M 

.071 

B 

.071 

100 

S 

.079 

M 

.083 

B 

.179 

110 

S 

2.96 

TIDAL 

M 

2.27 

CYCLE 

B 

3.40 

1 

120 

S 

5.60 

M 

5-97 

B 

6.4  5 

130 

S 

.583 

M 

• 535 

B 

.463 

V 

140 

M 

.052 

70 

M 

.012 

80 

M 

.150 

90 

S 

.047 

M 

.054 

B 

.058 

100 

S 

.063 

M 

.131 

B 

.219 

110 

S 

1.10 

M 

.950 

B 

1.09 

120 

S 

3.24 

M 

3.20 

B 

3.24 

TIDAL 

CYCI£ 

130 

S 

3.22 

3 

M 

3.26 

B 

3.10 

l4o 

S 

1.85 

M 

1.64 

B 

1.61 

150 

S 

• 553 

M 

.435 

B 

.307 

160 

M 

.148 

170 

M 

• 135 
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TEST  6 


APPENDIX  HI 


Correction  of  Observed  Concentrations. 

Tables  of  Corrected  Concentration. 

The  high  adsorption  of  the.  methylene  blue  on  any  surface  with 
which  it  came  in  contact  •==  the  plaster  and  concrete  of  the  model,  dust, 
sample  jars,  and  even  the  photometer  cells  ■==  caused  considerable  diffi- 
culty. It  was  found  in  laboratory  tests  that  the  adsorption  on  the  model 
itself  could  be  materially  reduced  by  painting  it  with  a clear  plastic 
paint  and  this  was  done  between  Test  1-A  and  the  rest  of  the  tests.  The 
residual  loss  on  the  model,  toge the r with  the  losses  on  other  surfaces,  was 
compensated  by  a correction  factor  calculated  in  the  following  way. 

If  no  dye  were  lost  then  the  integral  of  the  dye  concentration,  expressed 
as  weight  per  unit  volume,  over  the  entire  volume  of  the  model  would  equal 
the  weight  of  the  dye  introduced.  The  difference  between  this  integral 
and  the  original  weight  of  dye  is  therefore  a measure  of  the  dye  lost  by 
adsorption.  The  ratio  of  the  integral  to  the  weight  of  dye  introduced  may 
serve  as  a correction  factor  for  changing  observed  concentration  to 
corrected  concentration  which  is  an  estimate  of  the  value  the  observed 
concentration  would  have  shown  had  adsorption  not  been  present. 

Consideration  of  the  stations  where  samples  were  taken  across  the 
estuary  showed  that  after  the  first  one  or  two  tidal  cycles  the  mean 
sectional  concentration  could  be  estimated  from  the  mid~channel  samples. 
Using  the  mid-channel  observations,  concentration  was  plotted  against 
longitudinal  distance  and  a smooth  curve  drawn  through  the  points.  An 
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example  of  such  a curve  is  shown  in  Figure  A-1II-  1 where  the  open  circles 
denote  concentrations  estimated  from  single  samples  and  the  solid  circles 
denote  concentrations  estimated  from  an  average  of  samples  at  surface, 
middepth,  and  bottom.  In  drawing  the  smoothed  curves  more  weight  was 
given  to  these  latter  points  when  a departure  from  the  actual  estimates 
seemed  necessary. 

From  the  smoothed  curves  the  concentration  at  each  10,000  foot 
station  was  read.  Each  concentration  was  then  multiplied  by  the  appro- 
priate mean  cross  sectional  area  and  the  integration  was  carried  out 
numerically.  The  ratio  of  the  total  dye  represented  by  the  distribution 
curve  to  the  amount  of  dye  originally  introduced  was  formed.  This  ratio 
was  computed  for  each  tidal  cycle  and  plotted  as  a function  of  time.  A 
smooth  curve  through  these  points  provided  the  correction  factor. 

Four  correction  curves  of  this  sort  were  used,  one  for  Test  1-A, 
one  for  Tests  1 and  4,  one  for  Tests  Z and  3,  and  one  for  Tests  5,  6,  and  7. 

As  was  to  be  expected,  the  correction  factors  for  Test  I-A,  run 
before  the  model  had  been  painted,  were  significantly  higher  than  those 
for  Test  1,  run  after  the  model  had  been  painted,  although  the  tests  were 
otherwise  the  same. 

The  data  analyzed  in  the  text  of  this  report  have  all  been  corrected  in 
this  way.  When  the  difference  in  the  nature  of  the  decay  during  the  later 
tidal  periods  between  the  first  five  and  the  last  three  tests  was  noted  the 
question  arose  as  to  whether  the  difference  was  real  or  an  artifact  of  the 
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correction  method  Since  the  essential  difference  still  remains  if  raw 
data  are  used,  although  with  somewhat  different  values  for  the  logarithmic 
slopes,  it  seems  unlikely  that  the  difference  is  a result  of  the  applied 
correction.  Figure  A-III-2  shows  the  uncorrected  data  for  Tests  1-A 
and  6 as  an  example.  For  Test  1-A  the  logarithm  of  concentration  versus 
time  is  still  approximately  linear  while  for  Test  6 the  straight  line  is  still 
not  suitable.  Figure  A-III-2  should  be  compared  with  Figure  24  which 
shows  the  same  plot  for  the  corrected  data  from  the  same  two  tests. 

The  table  in  this  appendix  gives  the  corrected  concentrations  for 
each  tidal  cycle  by  10,000  foot  stations. 


i 


Station 


Corrected 
Cone . 


20 

0.05 

30 

0.17 

4o 

1.90 

TIDAL 

50 

10.30 

CYCLE 

53.5* 

10.56 

1 

6o 

6.09 

TO 

0.52 

8o 

0.03 

10 

0.05 

20 

0.20 

30 

0A7 

4o 

1.49 

TIDAL 

50 

5-57 

CYCLE 

57-5* 

6.83 

2 

6o 

6.77 

70 

3-31 

8o 

1.21 

90 

0.31 

10 

0.02 

20 

o.i4 

30 

0.38 

4o 

0.96 

50 

1.79 

TIDAL 

6o 

2.91 

CYCIE 

70 

4.24 

4 

75* 

4.48 

8o 

4.28 

90 

2.66 

100 

1.07 

no 

0.24 

20 

0.03 

30 

o.ll 

4o 

0.27 

50 

0.57 

TIDAL 

6o 

1.19 

CYCLE 

70 

2.05 

8 

80 

2.97 

90 

3.34 

92.5* 

3.35 

* Location  peak  concentration 
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TEST  3 


TEST  k 


| Station 

Corrected 

Cone . 

19C 

200 

TIDAL 

210 

CYCLE 

220 

0.14 

30 

230 

0.12 

240 

0.11 

(cont’d) 

250 

0.08 

0 

0.10 

10 

0.12 

20 

0.15 

30 

0.18 

4o 

0.25 

50 

0.36 

60 

0.50 

70 

O.56 

80 

O.58 

90 

0.73 

100 

0.86 

110 

0.99 

120 

1.04 

126* 

1.06 

130 

1.04 

TIDAL 

l4o 

0.94 

CYCLE 

150 

0.78 

4o 

160 

0.63 

170 

0.50 

180 

o.4i 

190 

0.33 

200 

0.30 

210 

0.2  6 

220 

0.25 

230 

0.23 

240 

0.18 

250 

0.17 

260 

0.15 

270 

0.13 

280 

0.12 

290 

0.10 

300  I 0.08 
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TEST  5 


TIDAL 

CYCLE 

30 


Station  Corrected 

Cone. 


220 

0.62 

2J0* 

0.64 

240 

0.57 

250 

0.52 

260 

0.46 

270 

0.43 

280 

0.38 

290 

0.33 

300 

0.28 

310 

0.16 

320 

0.08 

(cont'd) 
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TEST  7 


APPENDIX  IV 


A.  Auxiliary  Data  on  Salinity  Distribution,  High  and  Low  Water  Profiles, 
and  Fresh  Water  Inflow. 

B.  Table  of  Cross  Sectional  Area,  Inter se ctional  Volume,  and  Accumulated 
Volume  for  each  5000  ft.  Station. 


DELAWARE  RIVER  MODEL 
Salinity  Profile  Data 
Tests  2-3-U 
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